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C " 3 Abstract. We study the non homogeneous quadratic Bessel zeta function 

, C-Rs( s i v i a ) defined as the sum of the square of the positive zeros of the Bessel 

function J viz) plus a positive constant. In particular, we give explicit formulas 
for the main associated zeta invariants, namely poles and residua, Crb (0, v, a) 

and Cfl B (°. ^ a )- 
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From the point of view of differential geometry and mathematical physics, the 
Riemann zeta function appears as the operator zeta function associated to the 
Laplacian operator on the line segment ^Hj [HI © [IH ■ A natural generalisation 
of this setting, is to consider a Sturm Liouville operator instead, i.e. a singularity 
£f) [ at one of the end points |§] EH U [7] [HI El- This leads again to a concrete 

zeta function, namely the Bessel zeta function, where the sum is extended on the 
<-j | positive zeros of the Bessel function J„(z), and reduces for the opportune choice of 

O v to the classical Riemann case. Such a function was first considered and studied by 

Stolarsky in |21j , where formulas for poles and residua are given, and more recently 
by other authors, who calculated the associated zeta invariants by different methods 
PP |T2]. I n these notes, we study the non homogeneous version of this function. We 
determinate his poles and give formulas for the residua. In particular, we introduce 
two simple but quite general methods to calculate the value of the derivative at the 
origin, and therefore the regularized determinant of the associated Sturm-Liouville 
singular operator 0] [E] Q]I|. 
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Consider the constant singular Sturm Liouville operator 

2 d 2 4v 2 - 1 2 

on the line interval (0, 1], with positive real v and q (the null cases can be easily 
obtained as limit cases). L v has the discrete resolution [T3] 

.2 2 3v,n 2 i { \ y/^xJu 1 



I 2 ' Uv+l{jv,n) j 

where j„ ; „ are the positive zeros of the Bessel function J v {z) ordered in increasing 
way [22]. 
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For analogy with the Ricmann case, we consider the following non homogeneous 
quadratic Bessel zeta function defined by: 

oo 

C,rb{s,v,o) = ^2(nl >n + a 2 y%, 

n=l 

for Re(s) > 1, where nn^n — j VlTl and a is real and positive, and we study its 
analytical extension. We can state our main result for the Bessel zeta function, 
that comes as a corollary of the more general result state in Proposition El below 
concerning the zeta function associated to the operator L v + q 2 . 

Proposition 1. The function Crb(s, v, a) has an analytic extension to the complex 
s-plane, smooth in v, up to a discrete set of simple poles at s — 1, —1, —3, . . . , whose 
residua can be computed using the known asymptotic expansions for the Bessel 
functions ( more precisely, they are given by the residua of the function z(s, v, a, 7r) 
in Proposition^ multiplied by 2); in particular 

Resi {C, RB {s,v,a),s = 1)) = 1, 

Resi (( RB (s,v,a),s = -1)) = ~2~2 ~ \ ~ 71-2 ° 2 
This extension is regular at s = 0, and 

Res (Cfls(s, v, a), s = 0) = ~ (v + ~ 
Res (£ RB (s,v,a),s = 0) = ~ log Vz-r 1 ^™^ . 



Notice that the values for the homogeneous case, follow immediately using the 
series expansion for the Bessel function I v {z) for small z, namely (see \TG\) 

We will study the more general setting, i.e. the function 

oo 

z(s, q, I) = C(s; L v + q 2 ) = J^Kn + <tT* '■ 

n=l 

Convergence of the series for Re(s) > h follows from classical estimates on the 
zeros of Bessel functions 22 . and Ci?s(2s, v, a) = z(s, v, a, n). 

Wc present here two approaches to calculate the zeta invariants associated to 
z(s, v, q, I). In the first, we produce an analytic representation that can be effectively 
used to get all the invariants; this will be particularly useful to generalise the method 
for the calculation of the zeta invariants associated to a disk or to a cone j20]- In the 
second approach, we give a very general lemma (Lemma0J to deal with regularized 
products, and we apply it to calculate z(0, v, q, I) and z'(0, v, q, I). 
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Let start with the first approach. Using Mellin transform JSj [13 ■ we get the 
analytic representation 

1 f°° 

z(s, v, q, I) = — -r / t s f(t, v, q, l)dt, 



m Jo 

where the trace of the heat operator is 

oo 

f (t,v,q, I) =Tie- tL » = £V 

n=l 

and from this the complex representation 

f{t,u,q,l) = ^ [ e- M R(\,v,q,l)d\, 
2tt« Ja c 

where the contour is A c = {A € C | | arg(A — c)| = 7r/4}, for some < c < q 2 , and 
the trace of the resolvent is 



i?(A '^' = SA-(A^ + (j2 )- 



We now observe that it is easy to express such function in terms of special func- 
tions. In fact, taking logarithmic derivative of the infinite product representation 
of the Bessel function I v {z), we get |22| : 

Lemma 1. 

R{ ^,q,l) = ^-^ z \o g I v {lz), 



Here, z = \J q 2 — A, we set arg(g 2 — A) = on the line (— oo, q 2 ) and we fix the 
sector s_|_ = {z eC | | argz| < ir/2} for z. 

At this point it is worth observing that information about poles and residua 
of z(s,v,q,l) can be obtained using the representation introduced, asymptotics 
expansions for Bessel functions and classical arguments |15| . This is an easy 
way for producing the results relative to the so called 'constant case' when studying 
regular singular operators |H] [S]. More precisely, and for completeness, we can 
state the following results: 

Lemma 2. For small t, 

i 

f(t, v , q, i) = <*i(y, q, 0* (i " 1)/2 + o(i //2 ) = 

2 = 



-j: 



-1/2 



2 

where 



_ (-!)« lq 2i ^ ^ qV + 
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Proposition 2. The junction z(s,v,q,l) has an analytic extension to the whole 
complex s-plane, smooth in v and q, up to a discrete set of simple poles at s — 



i, — 5, — §,■■•, with residua 



r> / / n 1 A a 2 k(v,q,l) , n 1 „ 
Resi z(s, v,q, 1), s = — — k I = — r-j -r-, fc = 0, 1, 2, . 



2 ; r(i-fc) 

TTws extension is regular at s = 0, and 



Res o (z(s,z/,<z,0,s = 0) = -- \ v + - 

Notice that the contribution of the non homogeneity term q 2 is equally shared 
among all the poles; in other words, the homogeneous Bessel zeta function has the 
same poles, but with (possibly) different residua (see |21|). 

The information available is not enough to deal with the derivative, that is an 
harder point; though, we introduce the following quite more general purpose result. 

Lemma 3. Suppose the zeta function z(s,x) = XmLi a «( x ) ~ S has the following 
representation (everything smooth in x): 

z(s,x) = —^ t 8 - 1 — / e- xt R(X,x)d\dt, 
r(s) Jo 2m J Ac 

where the contour is as above, and the there is a primitive function — T(A, x) for 
the function R{\,x) = —4rT{\,x), satisfying the following properties: 

(a) T is analytic near A = 0. 

(b) for large A and fixed x, there is an asymptotic expansion 

T(A, i) = -.. + A(x) log(-A) + B(x) + .... 
Then, z(s,x) can be analytically extended at s = and 

z(0,x) = -A(x), 
z'(0,x) = -B(x)+T(0,x). 



Proof Integrating by part, first in A and hence in t, the given complex representation 
for z(s, x), we get 

s f°° If e~ At 
z(s,x) = — t*- 1 — —T(X,x)d\dt. 
T(s) J 2m J Ac -X 

At this point, it is convenient to isolate the constant part of T writing 

T(X, x) = Xi(A, x) + T (x); 

in fact, it is clear that Tq(x) gives no contribution, since 

-To{x)dX = 0. 



A C -A 

Next, since by definition the zeta function is well defined for large s, 
*(.,*) = ^ jVJL f ^T l{ X,x)dXdt + sif { s), 
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where / is regular near s = 0. Because of the pole at A = 0, we have to split the 
complex integral as follows to use the expansion for large A (small t): 



I = f -1=1 +Ti(0,x), 



where C c is a circle around the origin of radius c. Moreover, by assumption (b), 

f e _At f c~ A 

/ — -T 1 (X,x)dXdt = / —T 1 (X/t,v,q)dXdt = 

= • • • + jA{x) + A(x) ln(-A) + B(x) + . . . , 
where we have explicit-ed only the relevant part. This means that we can write 



z(s, x) = 



jA(x) - B(x) - -A(x) + Ti(0, x) 
s 



T(s + 1) 

where again g is regular near s — 0, and, from that, the thesis follows at once. □ 

Applying this argument to the function z(s, v, q, I), we have 

T(A, v,q,l) = u log Iz - log I v (Iz) - log 2T(i/ + 1), 
and hence we compute 



My, qJ) 



l 



i 



B(u,q,l) = -log27r + 



1 



iog?-i og 2 1/ r( ; / + i), 



T(0, v, q,l) = v log lq - log I„ - log 2T(i/ + 1). 



This gives: 



Rcs (z(s,i/,g,0,s = 0) = -- z/+ - , 



1 



and proves the following: 
Proposition 3. 



Reso(V(s, v, q, I), s = 0) = — log V 2-kI - . 



Let's turn to the second approach. First, we have the following: 

Lemma 4. Suppose two sequences a n and b n , n = 1,2,3,..., of real positive 
numbers are given, and satisfy the following conditions: 

(A) there are real s a , s& such that the two series 

oo oo 
n— 1 n— 1 

converge for Re(s) > s a , s& respectively; 

(B) ^/ie zeia function Cb(s) has an analytic extension at s = 

(C) |o n — 6„| < X definitely for some constant K, 

(D) s 6 < 1. 
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Then, the zeta function has an analytic extension at s = with £ a (0) 

Cb(0), the infinite product 



oc 



converges absolutely, and 



n—l 



c(o) = a(o)+io g c. 



Proof Let c„ = - 1. Then, |c„| < irt" 1 by (B) and (C) and hence £^° =1 c„ 
converges absolutely and so does the infinite product. We can assume c n < 1, and 
hence 



oo 

-S 



c-(«) = &(«)+£( fe )E & 

*!=i 

Now, the series 



s k 
n u n " 



V b~ s c k 



converges in a neighbourhood of s = for each k, by conditions (C) and (D). This 
means that £ Q can be analitycally extended to s — by using the extension of £b. 
In particular, evaluating at s = we get Ca(0) = Cf>(0) an d 



c(o) = a(o)+Ex T E & ^ 

fe=l ^ ' n=l 



5=0 fc=l 



E( j 8 



s=0 



fc-1 



= CUo) + E 5] c " = Cb(0) + E + c ») = Cb(o) + logC. 

□ 

We can apply this lemma to the present case as follows. Let a n = X 2 , n + q 2 , 

bn = jt \n + | (y — \Y\ 2 'i then, all the assumptions of Lemma0]are satisfied, and 
Cfe(s) is the Hurwitz zeta function (2s, \ (v + |)) . Thus, 

^(0) = 2io g rfi (v + l))-]o e 2K, 



2 , 

and we can compute (where u= ^{y — 5)): 

1 Ilri=l f 1 + 52^2 l+qip) ) ( 1 

— = lim 7 ; f- — lim 7 

c 2 ^o+ n°° hi ' 2 1 ^0+ n°° m 1 l2q2 \ m 1 . ' 2 



lim 



^(^)|r(«+i+ii)r ^ {lqY 



*-0+ I„(lq)T 2 (u + l) 7r"J„(l«)r 2 (| (!/+§)) 
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